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Abstract 



The process of a light neutral vector meson electroproduction is studied in the 
framework of QCD factorization in which the amplitude factorizes in a convolution 
fJL^ \ of the nonperturbative meson distribution amplitude and the generalized parton 

\q ■ densities with the perturbatively calculable hard-scattering amplitudes. We derive 

I . a complete set of hard-scattering amplitudes at next-to-leading order (NLO) for the 

| production of vector mesons, V = p°,uj,(f). 

> : 

(N . 
oo ■ 

1. The process of elastic neutral vector meson electroproduction on a nucleon, 

rn ■ 7*(g) N(p) -> V(q') N(p') , where V = p°, u, </>, (1) 

O ; 

& ' was studied in fix target and in HERA collider experiments. The primary motivation 

for the strong interest in this process (and in the similar process of heavy quarkonium 
production) is that it can potentially serve to constrain the gluon density in a nucleon 
[1, 2]. On the theoretical side, the large negative virtuality of the photon, q 2 = —Q 2 , 
j> | provides a hard scale for the process which justifies the application of QCD factorization 

methods that allow to separate the contributions to the amplitude coming from different 
scales. The factorization theorem [3] states that in a scaling limit, Q 2 — ► oo and xsj = 
Q 2 /2(p ■ q) fixed, a vector meson is produced in the longitudinally polarized state by 
the longitudinally polarized photon and that the amplitude of the process (1) is given 
by a convolution of the nonperturbative meson distribution amplitude (DA) and the 
generalized parton densities (GPDs) with the perturbatively calculable hard-scattering 
amplitudes. In this contribution we present the results of our calculation of the hard- 
scattering amplitudes at NLO. 

2. p 2 = p' 2 = m 2 N and q' 2 = m 2 M , where and tum are a proton mass and a meson 
mass respectively. The invariant cm. energy s = (q + p) 2 = W 2 . We define 

A = j/-p, P = V Afi t = A 2 , 

Q 2 

(q - A) 2 = m 2 M , x Bj = W2 + Q2 ■ (2) 

We introduce two light-cone vectors 

n 2 + = n 2 _ = , n + n_ = 1 . (3) 
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Any vector a is decomposed as 

a M = a + n+ + a~n1 + a±_ , a 2 = 2 a + a~ 
We choose the light-cone vectors in such a way that 



(4) 



p=(l + S)Wn+ + 



m 



N 



2(1 + £)W 
<m 2 N + A 2 



71- 



p' = (1 - OWn+ + 2(7-'^ U - + A± • (5) 
We are interested in the kinematic region where the invariant transferred momentum, 

/ 4 £2 



t = - 



2 m N 



(6) 



is small, much smaller than Q 2 . In the scaling limit variable £ which parametrizes the 
plus component of the momentum transfer equals £ = Xsj/(2 — %Bj)- 

GPDs are defined as the matrix element of the light-cone quark and gluon operators: 



2(Pn_) 



1 

" 2(Pn_) 



2mjv 



■«(p) 



(7) 



(Pn_) 7 2tt 



d^ e i X (Pz) n 



. a n_p(p>\G^ (-£) Gj (0 b)| 



z=An_ 



ft* (a, f , u(p') ^_ M (p) + f , t) u(p') 



, ia al3 n- a Ap 



2rriN 



u(p) 



(8) 



In both cases the insertion of the path-ordered gauge factor between the field operators is 
implied. Momentum fraction x, —1 < x < 1, parametrizes parton momenta with respect 
to the symmetric momentum P = (p+p')/2. In the forward limit, p' = p, the contributions 
proportional to the functions S q (x, £, t) and £ 9 (x, £, t) vanish, the distributions Ti q (x, £, t) 
and TC 9 (x, £, t) reduce to the ordinary quark and gluon densities: 



n q (x,0,0) = q(x) for x > , 
Ti g (x, 0, 0) = -q(-x) for a: < ; 
Ti 9 (x,0,0) = xg(x) for x > . 

Note that gluon GPD is even function of x, H 9 (x,^,t) = H 9 (—x,£,t). 
The meson DA <f>v(z) is defined by the following relation 



(9) 



1 

(Q\q(yhM-y)\ v L(p)) y ^o = p^fv Jdz< 



,i(2z-l)(py) 



<j>v(z) . 



(10) 



It is normalized to unity / (py(z)dz = 1 . Here z is a light-cone fraction of a quark, f v is 

o 

a meson dimensional coupling constant known from V — > e + e~ decay, f p = 198 ± 7MeV. 
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Al 7 2Ar- > y i iv = — — 7 dz<^v{z) y 



g y x) f»(x, e, t) + t (+) (z, (x, t)) + £ x ) F9(+) Z, t) 



(11) 



Here the dependence of the GPDs, DA and the hard-scattering amplitudes on factoriza- 
tion scale p F is suppressed for shortness. Since we consider the leading helicity non-flip 
amplitude, in eq. (11) the hard-scattering amplitudes do not depend on t. The account 
of this dependence would lead to the power suppressed, ~ t/Q, contribution, a is a fine 
structure constant, N c = 3 is the number of QCD colors. Qy depends on the meson flavor 
content. If one assumes it is -^(\uu) — \dd)), -^(\uu) + \dd)) and \ss) for p, uj and (j) 

respectively, than Q p = Q M = and Q$ — the sum in the last term of (11) runs 
over q = u, d for V = p, uj and q = s for V = and 

2 1 —1 
P p — p u — p p — _p^ — P 4> — t 

6u 3^2' d d 3V2 s 3 ■ 

F q ( + \x,£,t) = F q (x,£,t) - F q (-x,£,t) denotes a singlet quark GPD, F^\x,^,t) = 
J2 q =u,d,s F q ( + \x, £, t) stands for the sum of all light flavors. 

Due to odd C- parity of a vector meson <f>v(z) = 0y(l — z). Moreover, since V and 
7* have the same C- parities, 7* — > V transition selects the C-even gluon and singlet 
quark contributions, whereas the C-odd quark combination F q( ^~\x,^,t) = F q (x,£,t) + 
F q (—x,C,,t) decouples in (11). 

3. Below we present the results of our calculation of the hard-scattering amplitudes in 
the MS scheme. 

T q (z,x) may be obtained by the following substitution from the known results for a 
pion EM formfactor, see also [4], 

T q (z, x) = It (z, ^~ie)-T U ~ ie) } + {* - z} , (12) 



Avu \ 



<*s(/4) 

47T 



ci ( 2 [3 + ln(uu)] In 



'Q 



2 , + ln 2 (vu) 
Pp ) 



ln(v) ln(u) 28 \ n /5 , . , , (Q 2 
+61n(W) - - - — + A) o - 1 < VU ) ~ ln ( TT 



u 



V 

(yy^ -)- un?} 

+c 2 ( 2^- [Li 2 (u) - Li 2 (v) + Li 2 (v) - Li 2 (u) + \n(v) ln(w) - ln(w) ln(u)] 



v (f — -u) 3 

+2 [Li 2 («) + ^(u) - Li 2 (u) - Li 2 {v) + ln(u) \n(u) + ln(w) In(u)] 

2f-u) liivu 



, vu \n(vu) (v + u 
+A- r + 2— 



(f — u)'- 



(v — u)' 
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4 ln(f ) ln(-u) — 



(13) 



Here and below we use a shorthand notation u = 1 — u for any light-cone fraction. is 
a renormalization scale for a strong coupling, /3 = — is the effective number 
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of quark flavors. Cf 



Li 2 (z) = - f f ln(l - t). Also we denote 



Cl = C F , c 2 = Cf - 



T( + )(z,x) starts from NLO 



here 



f 2y + l 



(87r)zz 



2£ 



(14) 



(15) 



■" + 1 ln 2 (y + l)W 



R( Z ,y) + y(y + l) + ( l+Z y H{ A +{z ^_ zh 



yln(-y) + (y + l)ln(y + l) 
2/(2/ +1) 



(16) 



y + z (y + ^) 2 j 

where we introduced two auxiliary functions 

R(z, y) = z ln(-y) + z ln(y + l) + z \n(z) + z \n(z) , (17) 

H(z, y) = Li 2 {y + 1) - Li 2 {-y) + Li 2 (z) - Li 2 {z) + ln(-y) ln(z) - ln(y + 1) \n(z) . (18) 
For the gluonic contribution we obtain 



T a (z,x) 



«s(/4K 



zz (x + £ — it) (x — £ + it) 



1 «s(/4) ^ f _ a; - £ 



4. "M*' 



2e 



+ ie 



where 



ci (y\n(-y) (y + 1) ln(y + 1) 



(19) 



(20) 



+ In 



/4 



y + 1 



2( Cl -c 2 )(y 2 + (y+l) 2 ) 
2/(2/+ 1) 



+ ci (- + 2zm(z) 



( y ln(-y)-(y + l)ln(y + l)) + ^ 



- + ln(^) + ln(-y) + ln(y + 1) 
+ ( Cl (y(y + 1) + (y + z) 2 ) - c 2 (2y + l)(y + z)) 



i?(z, y) ln(-y) - ln(y + 1) + ln(z) - ln(z) 



(y + z)< 



+ 



2(y + z) 



y(y + l) + (y + z)- 
(y + z) 3 
Cl (2y+l)#(z,y) 
2(// • :! 



H(z,y) 
- 2ci 



-(c, - c 2 )(ln(^) - 2) + (y+1) ' n(j,+ 1) ) 4- c l( l 4- 3.) h,(*) 



V 2/ + J- 2/ 
(3 Cl -4c 2 ) /yln 2 (-y) , (y + 1) ln 2 (y + 1) ' 



y + 1 



+ 



y 



+ CiZ ln 2 (z) 



+ (ln(-y) + ln(y + 1)) ( Cl (zln(z) - + 2c 2 ) | + {z -> z} . 



4. Above equations give a complete description of a neutral vector meson electroproduc- 
tion with NLO accuracy. At leading order we reproduce known result [5], our results for 
the NLO correction are new. 

At high energies, W 2 ^> Q 2 , the imaginary part of the amplitude dominates. Leading 
contribution to the NLO correction comes from the integration region £ <C |x| <C 1, 
simplifying the gluon hard-scattering amplitude in this limit we obtain the estimate 



M--y* L N^V L N 

1 



/ 





dz 4>v( z ) 



zz 



2 i 7r 2 v^Lra asfvQv 

Q 2 zz 



N c Qi 



(21) 



F*(£,£,f) + — h 



J x 



IT 



Given the behavior of the gluon GPD at small x, F 9 (x,C,,t) ~ const, we see that NLO 
correction is parametrically large and negative unless one chooses the value of the fac- 
torization scale sufficiently lower than the kinematic scale. For the asymptotic form of 
meson DA, 4>y(z) = 6zz, the last term in (21) changes the sign at fiF = 7, for the DA 
with a more broad shape this happens at even lower values of //p. 
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